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Abstract 

Using CLEO data collected from CESR e+e~ collisions at the 4'{2S) resonance and nearby 
contimium at ^/s=S.67 GeV, we report the first significantly non-zero measurements of light vector- 
pseudoscalar hadron pair production (including pvr, lutt, pr], and K*^K^) and the 7T~^Tr~Tr^ final 
state, both from il^{2S) decays and direct e+e~ annihilation. 



*0n leave of absence from University of Chicago. 
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The fm puzzle poses one of the most enduring questions in strong interaction physics: 
why is the branching fraction for 4'{2S) — > pn at least twenty ^ times smaller than expected 
from scaling the J/ip pir rate by the ratio of dilepton branching fractions? The "12% 
rule" , a scaling conjecture generalizing this question for any decay mode, has as its underlying 
assumption that since charmonium decay to light hadrons must proceed through annihilation 
of the constituent cc into a photon or three gluons, the decay width should be proportional 
to the square of the cc wave function overlap at the origin. The rule's figure of merit is 

B{^{2S) ^ X)/B{J/^ ^ X) 

where B denotes a branching fraction and X a particular final state. Decays to dileptons also 
proceed via cc annihilation, and their branching fractions are well- measured jH , so their ratio 
makes a suitable denominator in Eq. (1). Several channels have Qx ~ 1 0, although some 
deviations from unity are expected The pn mode is not alone: significant suppressions 
also exist for at least one other vector-pseudoscalar (VP) channel {K*(892)~^K~) and three 
vector-tensor channels (pa2(1320), K*{S92)K*{U30), and o;/2(1270)) [l|,|3]. 

The continuing struggle to understand the pattern of Q-values has provoked many the- 
oretical explanations. For isospin-violating (IV) modes such as utt and pi], three-gluon 
mediated decay is suppressed, allowing the electromagnetic process of annihilation into a 
virtual photon to dominate. Whether Qjy ~ 1 remains a crucial open question. A recent 
review [2] of relevant theory and experiment concludes that none of the proffered theoret- 
ical explanations is satisfactory and also finds the underpinnings of the 12% rule overly 
simplistic. 

A major impediment to addressing the puzzle in a systematic manner is the dearth of 
ip{2S) branching fraction measurements. Experimental progress on key VP final states has 
remained dormant for many years. Continuum production, e"'"e~ ^ 7* — X, which is of 
interest in its own right is expected at levels that may affect 'ip{2S) backgrounds and 

will interfere ;6| with 4'{2S) decay, but has not yet been measured. Using e"'"e~ collision data 
acquired with the CLEO detector operating at the Cornell Electron Storage Ring (CESR), 
this Letter presents ip{2S) branching fractions and continuum cross sections for tt+tt^tt^; pvr, 
un, 4>n, pr], ur], 0r/, K*°{892)K0, K*+{892)K~] 6i(1235)7r. Where applicable, the inclusion 
of charge conjugate states is implied. We use p nn, vr" 77, u tt+tt'tt^, cj) K~^K~ , 
r/ ^ 77 and 7r+7r-7r°, K*^ K-7r+, K*+ K'^s^+ and K+tt^, and K^g 



TT ' 7T 



The CLEO 111 detector 0] features a solid angle coverage of 93% for charged and neutral 
particles. For the data presented here, the charged particle tracking system operates in a 
1.0 T magnetic field along the beam axis and achieves a momentum resolution of ~0.6% 
at p = 1 GeV/c. The cesium iodide (Csl) calorimeter attains photon energy resolutions of 
2.2% at i^^^ = 1 GeV and 5% at 100 MeV. Two particle identification systems, one based on 
ionization energy loss (dE/dx) in the drift chamber and the other a ring- imaging Cherenkov 
(RICH) detector, are used together to separate from vr^. The combined dE/dx-KlCH 
particle identification has efficiencies >90% and misdentification rates <5% for both tt"^ and 
K^. 

Half of the ip{2S) data and all the ^/s=?).Q^ GeV data were taken after a transition to 
CLEO-c in which CLEO Ill's silicon-strip vertex detector was replaced with a six- layer 
all-stereo drift chamber. The two detector configurations correspond to different accelera- 
tor lattices: the former with a single wiggler magnet and a center-of-mass energy spread 
AE=1.5 MeV, the latter (CESR-c |8|]) with the first half of its full complement (12) of 
wiggler magnets and Aii^=2.3 MeV. 
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The integrated luminosity (£) of the datasets was measured using e"'"e~ 77 
events j^. Event counts were normahzed with a Monte Carlo (MC) simulation based on the 
Babayaga event generator combined with GEANT-based detector modeling. The 
datasets have £=5.63 pb"^ on the peak of the ij{2S) (2.74 pb'^ for CLEO III, 2.89 pb~^ 
for CLEO-c) and 20.46 pb^^ at ^/s=3.67 GeV (all CLEO-c). The scale factor applicable to 
continuum yields in order to normalize them to ijj{2S) data, / = 0.268 ± 0.004, includes a 
2.6% correction to the C ratio to scale it by jlHf; the error includes both the relative 
luminosity and form factor s-dependence uncertainties. We also correct each final state's / 
for small efficiency differences between the ip{2S) and continuum samples caused by detector 
configuration. 

We base our event selection on charged particles reconstructed in the tracking system and 
photon candidates in the Csl calorimeter. Energy and momentum conservation is required 
of the reconstructed hadrons, which have momenta pi and total energy i?vis- We demand 
0.98 < E^ijy/s < 1.015 and - \p2\\/{y/s/c) < 0.02 (for tt+tt-ttO, pi = p^o and = 
p^+ + Pt^-), which together suppress backgrounds with missing energy or incorrect mass 
assignments. The experimental resolutions are smaller than 1% in scaled energy and 2% in 
scaled momentum difference. In order to suppress hadronic transitions to J/ip, we reject 
events in which any of the following fall within 3.05-3.15 GeV: the invariant mass of the 
two highest momentum tracks; or the recoil mass from the lowest momentum single 7r°, 
TT^vr" pair, or ti~^71~ pair. Feeddown from n'^n^J/ip, J/ip — >■ into 7r"'"7r~7r°, p~^n~, or 

{K^n'^)K~ is additionally suppressed by requiring M{fj,'^fj,~) < 3.05 GeV for those channels. 

MC studies were used to determine invariant mass windows for intermediate particle 
decay products. To reduce contamination from c<j/2(1270) [3| and ci;/o(600) [3| in ^i^r, we 
exclude Mt^tt <1.5 GeV. Similarly, pi] candidates with low mass rjir^ states are avoided 
with M(r77r^)min > 1.4 GeV. For 7r° —> 77, r] —>■ 77, and —>■ 7i~^7i~ candidates we use 
kinematically constrained fits of the decay products to the parent masses. Fake 7r° and r) 
mesons are suppressed with lateral shower profile restrictions and by requiring that their 
decays to 77 not be too asymmetric. 

For 7r"'"7r~7r°, p'^n~ , and prj {(pi]) with rj 77, one of the two final state charged particles 
must be positively identified as a tt^ [K^), but neither can be positively identified as a 
(tt^). Charged kaons in K*K must be identified as such, and any tt"^ candidate must not be 
identified as K^. Charged particles must not be identified as electrons using criteria based 
on momentum, calorimeter energy deposition, and dE/dx. The softer charged particle in 
two-track modes must have p < 0.425 x y^/c to suppress potential background from /i"'"/i~7 
in which a fake vr'' is found. Both tracks in two-track modes must satisfy | cos6'| < 0.83, 
where 9 is the polar angle with respect to the direction. 

The efficiency e for each final state is the average obtained from MC simulations of 



both detector configurations. The VP modes are generated [1^ with angular distribution 
(1 -|- cos^6') ^], 6i7r fiat in cos6', and n^n'Ti^ as in oj decay. We assume B{hi a;7r)=100%. 

Background contamination from other ip{2S) decays is determined from sidebands neigh- 
boring the signal windows in 7r°, rj, u, 0, Kg, K*, and bi candidate mass distributions. 
The sideband yields from the ip{2S) sample are decremented by the corresponding number 
of scaled continuum events (because scaled continuum events inside the signal window are 
subtracted separately) and by the small residual signal contributions expected, and then 
scaled to match the signal window size. 

We normalize the branching fractions to the total number of produced ip{2S) events. 
The technique described in Ref. [lH is applied to the datasets used here, resulting in a total 
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number of 4'{2S) decays of 3.08x10^. 

Kinematic distributions are shown in Figs. 1-4 and the event totals and efficiencies in 
Table I. We observe signals for several modes in both ip{2S) and continuum datasets. The 
significances S in the last column of Table I refiect the likelihood that the ip{2S) yields cannot 
be attributed to backgrounds alone. S is computed from trials in which Poisson fiuctuations 
of the ip{2S), continuum, and cross-feed contributions are all simulated to obtain a confidence 
level (CL) that a given mean ip{2S) signal /i combined with backgrounds would exceed or 
equal the observed event count. S is obtained from this procedure with /i=0. 

Table II shows the final results. We compute branching fractions with a straightforward 
subtraction of luminosity-scaled continuum yields; the value of the true branching fraction 
depends on the unknown three-gluon decay amplitudes and corresponding unknown phases. 
Statistical errors shown correspond to 68% CL and upper limits to 90% CL, and are obtained 
through simulated trials as described above. Values of Q are computed for each mode based 
on branching fractions from Ref. 0, except for B{J/ip tt+tt^tt^ ) = (2.10 ±0.12)% [13. 
Born-level cross sections at ^/s=3.67 GeV are also given and include an upward adjustment 
of 20% to account for initial state radiation fld\. 

The systematic errors on branching fractions share common contributions from the num- 
ber of produced 4'{2S) events (3%), uncertainty in / (1.5%), trigger efficiency (1%), electron 
veto (0.5% per veto), and MC statistics (2%). Other sources of uncertainty vary by channel; 
listed with their contribution to the systematic error, they stem from cross-feed subtractions 
(the change induced by ±50% cross-feed variation), accuracy of MC-generated polar angle 
and mass distributions (10% for bin, 14% for Tr+Tr^Tr*^), and imperfect modeling of charged 
particle tracking (1% per track), vr", rj and Kg finding (2% per vr*^ or rj, 5% per Kg), ir^/K^ 
identification (3% per identified ti/K), and mass resolutions (2%). Cross section system- 
atic errors include the above contributions, substituting an uncertainty in C (3%) for the 
normalization error, and accounting for uncertainties in the effects of initial and final state 
radiation (7%). Except for bin and h^h^tt^, statistical errors dominate. 

The ^'(25') results in Table II are consistent with previous measurements 0, where avail- 
able. Unlike other VP channels, the isospin-violating modes ujir and pr] are not strongly 
suppressed with respect to the 12% rule, an important new piece of the pn puzzle. The ratio 
B{ij{2S) K*+K-)/B{ij{2S) K*^K^) = O.Utofe is found to be much smaller than 
the equivalent ratio for J/ip decays, 1.19 ± 0.15 0|. Fig. 4 shows that ip{2S) tt+tt^tt^ 
decays have not only a distinct pvr component above the continuum contribution, but, unlike 
J/ip ^ TT^TT'n^ jl6||, which is dominated by pvr, also feature a much larger cluster of events 
near the center of the Dalitz plot. The pvr results reported here do not account for any 
cross-feed from this non-pvr component due to its uncertain source and shape. If five events 
inside the p mass window were attributed to the higher mass structure, the pir branching 
fraction would decrease by a quarter and its significance by one unit. 

The SU(3) expectation [17] for continuum cross sections is uott : prj : K*^K'^ : pre : (prj : 
K*+K- : cur/ : 07r = 1 : 2/3 : 4/9 : 1/3 : 4/27 : 1/9 : 2/27 : 0, in which a mixing angle 
9p satisfying sin6'p = —1/3, cos 6'p = 2-\/2/3 is chosen to describe rj-r]' mixing. With the 
striking exception of j^*Oj^o the measured VP continuum cross sections are consistent with 
Born-level calculations ja llSl and the above ratio predictions from SU(3). A least-squares 
fit for the common unit of cross section (corresponding to alun)), excluding K*^K^, yields 
(Tfit = 16.4 ± 2.7 pb with = 4.9 for 6 d.o.f.; a{K*^K^) exceeds (4/9)(Tfit by 3.0 standard 
deviations. Variations in 9p of ±10° induce changes of 15^ pb in o"fit- 

In summary, we have presented first evidence for ip{2S) decays to Tr^Tr^Tr*^, pix, prj, and 
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K*^K'^. Measurements for several other VP channels are also given. The results suggest 
that, for VP final states, ^'(25') decays through three gluons are severely suppressed with 
respect to the 12% rule and the corresponding electromagnetic processes are not. The decay 
ip{2S) — s> TT^TT^TT^ exhibits a p — > vrvr signal but has a much larger component at higher tttt 
mass. Continuum e+e~ cross sections for these final states are presented for the first time. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation and the U.S. Department of Energy. 

Note added in proof. Subsequent to the submission of this Letter, similar results from 
BES jlH became available. After correction for relative efficiencies and normalizations, the 
yields of events from ip{2S) and continuum datasets in the BES analyses are statistically 
consistent (within ilex) with those presented here. 
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TABLE I: For each mode: the efficiency, e; for ^/s=3.67 GeV data, the number of events, N^., 
and background from sidebands, TV^; for ^(2S) data, the number of events, N2S, the estimated 
continuum background, fN^, and background from other il){2S) decays, and the statistical 
significance S of the ip{2S) signal in units of a Gaussian standard deviation. 
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TABLE IL For each final state X: the branching fraction B{ip{2S) X), with statistical (68% 
CL) and systematic errors; the upper limit (90% CL), UL, for B; Q from Eq. (1); and a, the 
e+e~ — X Born-level cross section at ^/s=3.67 GeV. 
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FIG. 1: Distributions of scaled visible energy, E^\s/\/s, for labeled final states. Plots for pvr and ftivr 
sum over the charged and neutral states. Histogram entries are shown for ip{2S) data (points with 
error bars), scaled continuum (shaded histogram), and MC (dashed) with arbitrary normalization. 
The vertical arrows mark ends of signal selection ranges. 
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FIG. 2: Invariant mass distributions relevant to the final states indicated, one entry per event. 
Symbols are defined in Fig. 1. 
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FIG. 3: Invariant masses for the K*{892) (left) and Kg (right) candidates, one entry per event, 
for the final states indicated. Symbols are defined in Fig. 1. 
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FIG. 4: Distributions for the ir^n^ir^ final state: the Dahtz plots for (a) ip{2S) and (6) continuum 
data; (c) the vr+vr^, vr+vr'^, and tt^tt^ mass combinations (3 entries/event), and {d) the 77 mass. 
Symbols are defined in Fig. 1. 
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